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Abstract: We consider Ising-spin systems starting from an initial Gibbs measure v and evolving 
under a spin-flip dynamics towards a reversible Gibbs measure /j / v. Both v and /j, axe assumed 
to have a finite-range interaction. We study the Gibbsian character of the measure vS(i) at time 
t and show the following: 
' (1) For all v and /i, uS(t) is Gibbs for small i. 

(2) If both v and /j have a high or infinite temperature, then uS(t) is Gibbs for all t > 0. 
' (3) If 1/ has a low non-zero temperature and a zero magnetic field and fj, has a high or infinite 

temperature, then vS{t) is Gibbs for small t and non-Gibbs for large t. 

(4) If v has a low non-zero temperature and a non-zero magnetic field and [i has a high or infinite 
temperature, then vS(t) is Gibbs for small t, non-Gibbs for intermediate t, and Gibbs for large t. 
The regime where fj, has a low or zero temperature and t is not small remains open. This regime 
O I' presumably allows for many different scenarios. 



1 Introduction 



s 

> 

^ ■ Changing interaction parameters, like the temperature or the magnetic held, in a thermody- 

namical system is the preeminent way of studying such a system. In the theory of interacting 
particle systems, which are used as microscopic models for thermodynamic systems, one as- 
sociates with each such interaction parameter a class of stochastic evolutions, like Glauber 
dynamics or Kawasaki dynamics. 

In recent years there has been extensive interest in the quenching regime, in which one 
starts from a high- or inhnite-temperature Gibbs state and considers the behavior of the 
system under a low- or zero-temperature dynamics. This is interpreted as a fast cooling 
procedure (which is different from the slow cooling procedure of simulated annealing) . One is 
interested in the asymptotic behavior of the system, in particular, the occurrence of trapping 



in metastable frozen or semi- frozen states (see |ftl |, 34], [|35f[ , |i2| , 33|1 , [36| , 
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Another regime that has been intensively studied is the one where, starting from a low- 
non-zero-temperature Gibbs state of Ising spins in a positive magnetic field, one considers a 
low-non-zero-temperature negative-magnetic-field Glauber dynamics (see [^] and references 
therein). Under an appropriate rescaling of the time and the magnetic-field strength, one 
finds a metastable transition from the initial plus-state to the final minus-state. 

In this paper we concentrate on the opposite case of the unquenching regime, in which one 
starts from a low-non-zero-temperature Gibbs state of Ising spins and considers the behavior 
of the system under a high- or infinite-temperature Glauber dynamics. This is interpreted as 
a fast heating procedure. As far as we know, this regime has not been studied much (see e.g. 
Up), as no singular behavior was expected to occur. Although we indeed know that there is 
exponentially fast convergence (cf. |p3fl . Chapter 1, Theorem 4.1, and |}^]) to the high- 
or infinite-temperature Gibbs state (i.e., the asymptotic behavior is unproblematic) , we will 
show that at sharp finite times there can be transitions between regimes where the evolved state 
is Gibbsian and regimes where the evolved state is non-Gibbsian. 

In the light of the results in |J, Chapter 4, on renormalization-group transformations, it 
should perhaps not come as a surprise that such transitions can happen. Indeed, we can view 
the time-evolved measure as a kind of (single-site) renormalized Gibbs measure. Even though 
the image spin o~t(%) at time t at site x is not a (random) function of the original spins o~o(y) 
at time for y in only a finite block around x, by the Feller character of the Glauber dynamics 
it depends only weakly on the spins o~o(y) with y large. In that sense the time evolution is 
close to a standard renormalization-group transformation, without rescaling, and so we can 
expect Griffiths-Pearce pathologies. 

We will prove the following: 



(1) For an arbitrary initial Gibbs measure and an arbitrary Glauber dynamics, both having 
finite range, the measure stays Gibbs in a small time interval, whose length depends on 
both the initial measure and the dynamics (Theorem 4.1). This result, though somewhat 



surprising, essentially comes from the fact that for small times the set of sites where a 
spin flip has occurred consists of "small islands" that are far apart in a "sea" of sites 
where no spin flip has occurred. 

(2) For a high- or infinite-temperature initial Gibbs measure and a high- or infinite-tempera- 
ture Glauber dynamics, the measure is Gibbs for all t > (Theorems 5.11 and |6.15| ). 



(3) For a low-non-zero-temperature initial Gibbs measure and a high- or infinite-temperature 
Glauber dynamics, there is a transition from Gibbs to non-Gibbs (Theorems 5.16 and 
6.18). This result is somewhat counter-intuitive: after some time of heating the system 
it reaches a high temperature, where a priori we would expect the measure to be well- 
behaved because it should be exponentially close to a Completely Analytic (see J?]]) 
high-temperature Gibbs measure. As we will see, this intuition is wrong. However, from 
the results of |29] it follows that this transition does not occur when the initial measure 
is a rigid ground state (zero-temperature) measure (i.e., a Dirac measure). 

(4) For a low-non-zero-temperature initial Gibbs measure and a high- or infinite-temperature 
Glauber dynamics, there possibly is a transition back from non-Gibbs to Gibbs when 
the Hamiltonian of the initial Gibbs measure has a non-zero magnetic field (Theorems 
516] and gUp. 
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The complementary regimes, with a low- or zero-temperature Glauber dynamics acting over 
large times, are left open. 

In Section |2| we start by giving some basic notations and definitions, and formulating some 
general facts. 

In Section |3] we give representations of the conditional probabilities of the time-evolved 
measure and clarify the link between the Gibbsian character of the time-evolved measure and 
the Feller property of the backwards process. These results are useful for proving the "positive 
side", i.e., for showing that the time-evolved measure is Gibbsian. We use a criterion of |J, 
Chapter 4, Step 1, or [10] to identify bad configurations (points of essential discontinuity of 
every version of the conditional probabilities) as those configurations for which the constrained 
system (i.e., the measure at time conditioned on the future bad configuration at time t > 0) 
exhibits a phase transition. This criterion will serve for the "negative side", i.e., for showing 
that the time-evolved measure is non-Gibbsian. 

In Section ||] we prove that for an arbitrary initial measure and an arbitrary dynamics, 
both having finite-range interactions, the measure at time t is Gibbs for all t G [0, to]) where 
to depends on the interactions. 

In Section || we treat the case of infinite-temperature dynamics, i.e., a product of inde- 
pendent Markov chains. This example already exhibits all the transitions between Gibbs and 
non-Gibbs we are after. Moreover, it has the advantage of fitting exactly in the framework of 
the renormalization-group transformations: the time-evolved measure is nothing but a single- 
site Kadanoff transform of the original measure, where the parameter p(t) of this transform 
varies continuously from p(0) = oo to p(oo) = 0. For the case of a low-temperature initial 
measure we restrict ourselves to the cf-dimensional Ising model. 

In Section |6] we show that the results of Section [| also apply in the case of a high- 
temperature dynamics. The basic ingredient is a cluster expansion in space and time, as 
developed in [^] and worked out in detail in | ]25[j . This is formulated in Theorem 6.3 and is 
the technical tool needed to develop the "perturbation theory" around the infinite-temperature 
case. 

In Section fj] we give a dynamical interpretation of the transition from Gibbs to non-Gibbs 
in terms of a change in the most probable history of an improbable configuration. We show 
that the transition is not linked with a wrong behavior in the large deviations at fixed time, 
and we close by formulating a number of open problems. 



2 Notations and definitions 

2.1 Configuration space 

The configuration space of our system is Q = { — 1, +l} zd , endowed with the product topology. 
Elements of f2 are denoted by a, rj. A configuration a assigns to each lattice point x G Z d a 
spin value cr(x) G {— 1, +1}. The set of all finite subsets of Z d is denoted by S. For A G S 
and a G Q, we denote by a\ the restriction of a to A, while Q\ denotes the set of all such 
restrictions. A function / : — > R is called local if there exists a finite set A C Z d such that 
f{rj) = f(o~) for a and r\ coinciding on A. The minimal such A is called the dependence set 
of / and is denoted by Df, The vector space of all local functions is denoted by C. This is 
a uniformly dense subalgebra of the set of all continuous functions C(f2). A local function 
/ : — > R with dependence set Df C A can be viewed as a function on With a slight 
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abuse of notation we use / for both objects. For a, r\ e $7 and A C Z d , we denote by atj]K c the 
configuration whose restriction to A (resp. A c ) coincides with a a (resp. t/a c )- For x6Z d and 
a e Q, we denote by r x cr the shifted configuration defined by r x a(y) = a(y + x). A sequence of 
probability measures fi\ on £l\ is said to converge to a probability measure ^ on (notation 
MA -> /«) if 



2.2 Dynamics 



lim fdfJL A = fdfi V/ e £. (2.1) 



The dynamics we consider in this paper is governed by a collection of spin-flip rates c(x, a), 
x G Z d , uefi, satisfying the following conditions: 

1. Finite range: c x : a i— > c(x, cr) is a local function of <r for all x, with diam(.D Ca ,) < R < oo. 

2. Translation invariance: t x cq = c x for all x. 

3. Strict positivity : c(x, <r) > for all x and a. 

Note that these conditions imply that there exist e, M e (0, oo) such that 

< e < c(x, <r) < M < oo Vx e Z d , (2.2) 
Given the rates (c^), we consider the generator defined by 

Lf = CxVJ V/ e C, (2.3) 

x&lA 

where 

V x /((7) = /(O - /(a). (2.4) 

Here, a 3 ^ denotes the configuration defined by o~ x (x) = —a{x) and a x (y) = o~(y) for y ^ x. 
In [^3|, Theorem 3.9, it is proved that the closure of L on C(Q) is the generator of a unique 
Feller process {at : t > 0}. We denote by S(t) = exp(iL) the corresponding semigroup, by F a 
the path-space measure given cr = cr, and by E CT expectation over P^. 

A probability measure [i on the Borel <r-field of 0, is called invariant if 



L/d^ = V/ e C. (2.5) 



y (Lf)gdfi = J f(Lg)d(i Vf,g G C. (2.6) 



It is called reversible if 



Reversibility implies invariance. For spin-flip dynamics with generator L defined by (|27 
reversibility of [x is equivalent to 



c(x,a x )-^ = c(x,a) \/x£Z d ,a£tt, (2.7) 
d/j, 



where \x x denotes the distribution of <r x when a is distributed according to /x. Note that (|2T 
implies the existence of a continuous version of the Radon-Nikodym (RN)-derivative dfj, x jd\i. 
This will be important in the sequel. 
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2.3 Interactions and Gibbs measures 



A good interaction is a function 

f/:5xfl^R, (2.8) 
such that the following two conditions are satisfied: 

1. Local potentials in the interaction: U(A,o~) depends on a(x),x G A, only. 

2. Uniform summability: 

V sup | [/(A, a) | < oo x£Z d . (2.9) 

The set of all good interactions will be denoted by B. A good interaction is called translation 
invariant if 

U(A + x, T_ x a) = U(A, a) VA 6 S, x G Z d , a G fl (2.10) 

The set of all translation- invariant good interactions is denoted by Bu- An interaction U 
is called finite-range if there exists an R > such that U(A, a) = for all ^4 G 5 with 
diam(A) > i?. The set of all finite-range interactions is denoted by B^ r and the set of all 
translation-invariant finite-range interactions by B{[ . For U G B, £ G 12, A G S, we define the 
finite-volume Hamiltonian with boundary condition £ as 

4w= E ^(a^aCac) (2.H) 

AnA^0 

and the Hamiltonian with free boundary condition as 

4W = j;[/(A4 (2.12) 

AcA 

which depends only on the spins inside A. Corresponding to the Hamiltonian in ( |2.11| ) we 
have the finite-volume Gibbs measures /i» , A G S, defined on Q by 

/ /KJ/tfta) = £ /K<A.) eXP[ -J iW1 . (2-13) 

where denotes the partition function normalizing to a probability measure. 

For a probability measure /ionfi, we denote by f£ the conditional probability distribution 
of a(x),x G A, given a\c = £ A c. Of course, this object is only defined on a set of /U-measure 
one. For A G S,T G S and A C T, we denote by /ir(oA.|C) the conditional probability to 
find a\ inside A, given that £ occurs on T \ A. For U G B, we call /z a Gibbs measure with 
interaction C7 if its conditional probabilities coincide with the ones prescribed in ( |2.13| ), i.e., 
if 

/i£=/x^' C fi-a.s. AG5,CGfi. (2.14) 

We denote by Q{U) the set of all Gibbs measures with interaction U. For any U G £>, Q{U) is 
a non-empty compact convex set. The set of all Gibbs measures is 

G = |J £([/). (2.15) 
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Note that Q is not a convex set, since for U and V in Bu, convex combinations of Q{U) and 
GiV) are not in Q unless Q(U) = Q{V) (see ]|] section 4.5.1). 

Remark: We will often use the notation H = ^2 A U(A, •) for the "Hamiltonian" correspond- 
ing to the interaction U. This formal sum has to be interpreted as the collection of "energy 
differences", i.e., if a and rj agree outside a finite volume A, then: 

H{r,)-H{a)= £ [U(A, rj) - U(A, a)]. (2.16) 



Definition 2.17 A measure ji is called Gibbsian if ji £ Q, otherwise it is called non- 
Gibbsian. 



2.4 Gibbsian and non-Gibbsian measures 



In this paper we study the time-dependence of the Gibbsian property of a measure under the 
stochastic evolution S(t). In other words, we want to investigate whether or not uS(t) £ Q at 
a given time t > 0. 



Proposition 2.18 The following three statements are equivalent: 

1. n e Q. 

2. admits a continuous and strictly positive version of its conditional probabilities 
Ae5,(Gn. 

3. ji admits a continuous version of the RN- derivatives d\i x /dfi, x G Z d . 



Proof. See [2| and f}9|. 



We will mainly use item 3 and look for a continuous version of the RN-derivatives d[i x /dfj, by 
approximating them uniformly with local functions. 

A necessary and sufficient condition for not to be Gibbsian (fj, £ Q) is the existence of a 
bad configuration, i.e., a point of essential discontinuity. This is defined as follows: 

Definition 2.19 A configuration rj £ Q is called bad for a probability measure ji if there 
exists e > and x £ Z rf such that for all A £ S there exist T D A and £, £ £ Q, such that: 

\nr(.o-(x)\rj A \ {x} ( T \ A ) - Hr(o-(x)\r] A \ {x} ^ r \ A )\ > e. (2.20) 

Note that in this definition only the finite-dimensional distributions of \x enter. It is clear that 
a bad configuration is a point of discontinuity of every version of the conditional probabilities 
of ji. Conversely, a measure that has no bad configurations is Gibbsian (see e.g. [p?|l ). 



2.5 Main question 

Our starting points in this paper are the following ingredients: 

1. A translation invariant initial measure v £ G(U U ), corresponding to a finite-range 
translation-invariant interaction U u £ as introduced in Section 2.3. 



6 



2. A spin-flip dynamics, with flip rates as introduced in Section This dynamics has 
a reversible measure /u, which satisfies 

"-f - ^4- (2-21) 

d\x c(x,a x ) 



Hence, by Proposition 2.18 there exists an interaction f7„ G B such that fi G G(U^). Since 
the rates are translation invariant and have finite range, this interaction can actually be 
chosen in B{[ and satisfies (recall (|2,11 ) and ( |2.14| )) 



exp(^[^(A^)-^(A^)]J • (2.22) 



Without loss of generality we can take the rates c(x, a) of the form 



C (x,a)=e X pQ^[C/ M (A,a)-C7 M (A^)]j • (2.23) 



A finite-volume approximation of the rates in ( |2.23| ) that we will often use is given by 

c A (x, a) = exp [H£(a) - H»(a x )] , (2.24) 

where H 1 ^ is the Hamiltonian with free boundary condition associated with the interaction 
(recall ( |2,12j )). These rates generate a pure-jump process on fix = { — 1> with generator 

(L A /)(-) = ca(x, .)V«/(.) V/ G C. (2.25) 

Since L^f converges to Lf as A f 7j d for any local function / G C, the corresponding semigroup 
S\ (t) converges strongly in the uniform topology on C(Q) to the semigroup S(t), i.e., S\(t)f — ► 
S(t)f as A | Z d in the uniform topology for any / G C(O). Therefore we have the following 
useful approximation result. Let i/bea probability measure on Q and z^a its restriction to f^A 
(viewed as a subset of O). Then 

lim i/aSaCO = "5(t), (2.26) 



where the limit is in the sense of (2.1). If v G Q(U V ) is a Gibbs measure, then we can replace 



the finite-volume restriction i>a by the free-boundary-condition finite-volume Gibbs measure 
(in the case of no phase transition), or by the appropriate finite- volume Gibbs measure with 
generalized boundary condition that approximates v (in the case of a phase transition). 

The main question that we will address in this paper is the following: 



Question: 

Is vS(t) = Vf a Gibbs measure? 

In order to study this rather general question we have to distinguish between different regimes, 
as defined next. 
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Definition 2.27 U € B is a high-temperature interaction if 



sup J2(\A\ -1) sup \U(A,a) - U(A,a')\ < 2. (2.28) 



CT,cr'eSl 



Equation ( |2.28 ) implies the Dobrushin uniqueness condition for the associated conditional 



probabilities $ , A G 5,( £ O. In particular, it implies that = 1 (i.e., no phase 

transition). Note that it is independent of the "single-site part" of the interaction, i.e., of the 
interactions U({x},a). 

Remark: We interpret the above norm as an inverse temperature, so small norm means high 
temperature. 

Definition 2.29 We call: 



1. an initial measure v "high-temperature" if it has an interaction satisfying i\2.28 ), and 
write T v » 1. 

2. an initial measure v "infinite-temperature" if it is a product measure, (i.e., if the cor- 
responding interaction U u satisfies U u (A,a) = for all A with \A\ > 1), and write 
T v = oo. 

3. a dynamics "high-temperature" if the associated reversible Gibbs measure \x has an in- 



teraction Un satisfying ( 2.28 ), and write >> 1. 



4- a dynamics "infinite-temperature" if the associated reversible measure \x is a product 
measure (i.e., if the corresponding interaction satisfies U^(A,a) = for all A with 
\A\ > 1), and write = oo. 

As we will see in Section ||, the study of infinite-temperature dynamics is particularly instruc- 
tive, since it can be treated essentially completely and already contains all the interesting 
phenomena we are after. 



3 General facts 

3.1 Representation of the RN-derivative 



As summarized in Proposition |2.18| , an object of particular use in the investigation of the 
Gibbsian character of a measure is its RN-derivative d[i x jd[i w.r.t. a spin flip at site x. In 
this section we show how to exploit the reversibility of the dynamics in order to obtain a 
sequence of continuous functions converging to the RN-derivative of the time-evolved measure 
vi = vS{t) w.r.t. spin flip. Let us first consider the finite-volume case. We start from the 
finite-volume generator 



LaH<t) = CA(x,a)(f(a x ) - f(a)), (3.1) 



where the finite-volume rates c\(x, •) are given by ( |2.24 ). Suppose that our starting measure 



v £ G{U U ) is such that |C?(t7j,)| = 1, which implies that the free-boundary-condition finite- 
volume approximations v\ converge to v. The free-boundary-condition finite-volume Gibbs 



S 



measure fi A , corresponding to the interaction is the reversible measure of the generator 
L A . We can then compute, using reversibility, 



du A s A (ty 

du A S A (t) 



fdu A s A (tr^\ / dfi A s A {t) x ^\ Ufi A s A (t) ^ 



\d/j, A S A (t) 
( dv A S A {t) 
\dfi A S A (t) 



V dfi A S A (t) 



, , i , A d/j, A S A (t) 



d/x A / \dv A S A (t) 



\di^ A S A (t) 
(<0 



a>A 



(3.2) 



Definition 3.3 H^ v (a) = XmcaI^m^' ^ ~~ ^(Ao~)]- -Woie i/iai i/iis "difference Hamilto- 
nian" depends on both the initial measure and the dynamics. 



Using this definition, we may rewrite 



as 



dv A S A (t)* = dpi K* (expfO*)]) 



du A S A (t) 



<W ' EA (exp[^r («•*)]) 



(3.4) 



where E^ denotes the expectation for the process with semigroup S A (t) starting from a. Since 
this semigroup converges to the semigroup S(t) of the infinite- volume process as A — ► Z rf , we 
obtain the following: 



Proposition 3.5 For any a £ Q and t > 0, 

dvS(t) 



dvS{t) 



((T) dfx ATZ^ E CT (exp[H^(a t )}) 



(3.6) 



where this equality is to be interpreted as follows: if the limit in the RHS of \3.(fy is a limit 
in the uniform topology, then it defines a continuous version of the LHS. 



Proof. The claim follows from a combination of ( |2.26 ) and (3^) with Lemma 3/7 below. 



Lemma 3.7 If ' v n — > v weakly as n — » oo, and dv x /dv n £ C(Q) existe /or any n G N and 
converges uniformly to a continuous function ^> , then 



hm^ 

n]oo dv n 



dv x 
dv 



Proof. Let / : f2 — ► R be a continuous function. Define X : f2 
f o X : 0, — > R is a continuous function. Therefore 



(3.8) 

O by 0a;(o") = ct 3 ^. Then also 



{f o6 x {a))u(da) 
lim / (fo6 x (a)) v n (do) 
lim / -^(a)f(o-)Mdo-) 
lim / <&(a)f(a)v n (da) 

njoo J 



(3.9) 
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where the fourth equality follows from 

lim f _ * (<T ) f(a)u n (da) < lim ||/U|^ - *[U = 0. (3.10) 



Since ( |3.9| ) holds for any continuous function /, the statement of the lemma follows from the 
Riesz representation theorem. I 



Proposition [3^, combined with Proposition [2.18| , will be used in Sections ^-|6| to prove 
Gibbsianness. 



3.2 Path-space representation of the RN-derivative 

An alternative representation of the RN-derivative dvf jdvt is obtained by observing that 
vt = uS(t) is the restriction of the path-space measure Pj, to the "layer" {t} x O. In 
some sense, this path-space measure can be given a Gibbsian representation with the help 
of Girsanov's formula. The "relative energy for spin flip" of this path-space measure is a 
well-defined (though unbounded) random variable. Conditioning the path-space measure RN- 
derivative for a spin flip at site x in the layer {£} x tt, we get the RN-derivative dvf jdvt- More 
formally, let us denote by ir t the projection on time t in path space, i.e., 7r f (u;) = u>t with 
lo £ D([0,t],Q) the Skorokhod space. By a spin flip at site x in path space we mean a 
transformation 

Qx ■ D([0, t],fl) — > D([0, £], Q) (3.11) 

such that 

(7r t (u;)r = 7T t (@ x (uj)). (3.12) 
Different choices are possible, but in this section we choose 

rn> i \\i \ / -v(s,x) for y = x, < s < t, , . 

(e.( W ))( a> y) = | uM otherw . se _ (3.13) 

Let J-m denote the cr-field generated by the projection irt- Then we can write the following 
formula: 

This equation is useful because of the Gibbsian form of the RHS of fljTjj ) given by Girsanov's 
formula, as shown in the proof of the following: 

Proposition 3.15 Let v be a Gibbs measure on $7. For any t > 0, 

vS{t) x « uS(t) (3.16) 
and the RN-derivative can be written in the form 



duS{t) v 



dv x 

O 7T W x | T[ t ] 



(3.17) 



dv 

where ^f x : D([0,t],Q) — > R is a continuous function on path space (in the Skorokhod topology). 
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Proof. We first approximate our process by finite-volume pure-jump processes and use Gir- 
sanov's formula to obtain the densities of these processes w.r.t. the independent spin-flip 
process. Indeed, denote by P A the path-space measure of the finite-volume approximation 



with generator (|2,25|) and by P^' the path-space measure of the independent spin-flip process 
in A, i.e., the process with generator 



(3.18) 



We have for / : f2 



such that Dj C A 



/ 



f{a) vS{tf{da) 



lim / u(da) / F*(duj)f(7r t (@ x (u))) 



lim / v{da) I F^°(du;)^(u) f {7Tt { q x{uj))) . (3 . 19) 



d. 



dK 



Since P CT ' is the path-space measure of the independent spin-flip process, the transformed 
measure P^' ° & x equals P^i°. Abbreviate 



Fa(u) 



jA 

■ wo 



(3.20) 



Then we obtain 



v{da) I F£'°(cLj)F A (u)f(iT t (®x(u))) 



/jpA,0 
P^(^)^(a;)F A (e i .M)/(7r < H) 

du x f 



(3.21) 



where can be computed from Girsanov's formula (see |^4|] p. 314) and for A large enough 
reads 



^x,a( u ) = ex P 



[ t log^^\dN!(u)+ (\c(y^s)-c{y^ s )]ds 



\y— x\<R ' 



, (3.22) 



where Nf{ui) is the number of spin flips at site y up to time t along the trajectory uj. We thus 
obtain the representation of (|3.17| ) by observing that ^ x K does not depend on A for A large 
enough and using the convergence of P A to P CT as A f 7L d . Indeed, the only point to check is 
that 

' dv x \ 

— ott \ V x £L\F V ), (3.23) 

so that the conditional expectation in ( 3.17p is well-defined. However, this is a consequence 
of the following two observations: 



1. dv x jdv is uniformly bounded because v € Q. 
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2. For fy x we have the bound 

\* x (")\<e 2Ct (j) , (3-24) 

where, as before, M and e are the maximum and minimum rates, N t ' (lo) is the total 
number of spin flips in the region {y : \y — x\ < R} up to time t along the trajectory to. 
Since the rates are bounded from above, the expectation of the RHS of (|3.24| ) over 
is finite uniformly in a. 



3.3 Backwards process 



Proposition 3.15 provides us with a representation of the RN-derivative dvf jdvt that can be 
interpreted as the expectation of a continuous function on path space in the backwards process. 
The backwards process is the Markov process with a time-dependent transition operator given 
by 

(T v (a,t)f)(<r)=E v (fon a \a t = <r) 0<s<t, (3.25) 

where (-|ot = a) is conditional expectation with respect to the <r-field at time t. Notice that 
this transition operator depends on the initial Gibbs measure v and is a function of s and 
t (time-inhomogeneous process). Although the evolution has a reversible measure fj,, at any 
finite time the distribution at time t is not [i. This causes essential differences between the 
forward and the backwards process. 

The dependence of T u (s, t) on v is crucial and shows that even for innocent dynamics, like 
the independent spin-flip process, the transition operators of the backwards process may fail 
to be Feller for certain choices of v (see Section |B| below). In general, the independence of the 
Poisson clocks that govern where the spins are flipped (in the backwards process this means 
were flipped) is lost. 

In order to have continuity of the RN-derivative dvf /dist, it is sufficient that the operators 
T u (s,t) have the Feller property, i.e., map continuous functions to continuous functions. 

Proposition 3.26 // v is a Gibbs measure, then: 

T u {s, i)C(O) C C(Q) V < s < t < t => vS{t) G Q V <t <t . (3.27) 



Proof. This is an immediate consequence of Proposition 3.15. See also [20]. 



As in Section 3.1, we can thus hope to approximate the transition operators of the back- 



wards process by "local operators" (operators mapping £ onto £). 
Proposition 3.28 For any a G fl and < s < t, 

(T v (s,t)f){a) = hm y — , (3.29) 

ATZ d E CT (exp [H% (a t )}) 



where this equality is to be interpreted as follows: if the limit in the RHS of ( 3.29 ) is a limit 
in the uniform topology, then it defines a continuous version of the LHS. 
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Proof. Let us first compute T u (s,t) in the case of the finite- volume reversible Markov chain 
with generator ( 2. 25] ). For the sake of notational simplicity, we omit the indices A referring 
to the finite volume, and abbreviate v a = vS(s): 



(T v (s,t)f)((r) = ^2pts(v,° 



v t {a) 



«*> (I) w 



A*a(»7)' 

v 



S(s) 



fo) 



/fa) 



5(t 



S(t) 



E CT (exp[iT£<> t )]/(a t _ s )) 



E CT (exp[iZ£> t )]) 



(3.30) 



where H^' u is defined in Definition 



Propositions 3.26| and 3.2S| are the analogues of Propositions 2.18 and 3.5. We will not 
actually use them, but they provide useful insight. 



3.4 Criterion for Gibbsianness of uS(t) 

A useful tool to study whether i/S(t) 6 Q is to consider the joint distribution of (o-Q,a t ), 
where do is distributed according to v. Let us denote this joint distribution by Pt, which can 
be viewed as a distribution on {—1, +1} S with S = Z rf © Z d consisting of two "layers" of Z d . 
The correspondence between v t and uS{t) is made explicit by the formula 



MdcT,drj)f(a)g(rj) = / u(da)(fS(t)g)(a) f,g€C. 



(3.31) 



Now, for reasons that will become clear later, z?t has more chance of being Gibbsian than 
uS{t). The latter can then be viewed as the restriction of a Gibbs measure of a two-layer 
system to the second layer. Restrictions of Gibbs measures have been studied e.g. in |]37fl , 
HH Q, [|^], @, and it is well-known that they can fail to be Gibbsian, and most examples 
of non-Gibbsian measures can be viewed as restrictions of Gibbs measures. Formally, the 
Hamiltonian of v% is 

H t (a, n) = H v (a) - logger, r,), (3.32) 

where pt(cr,rj) is the transition kernel of the dynamics. Of course, the object log pt(cr, rj) has 
to be interpreted in the sense of the formal sums U {A, a) introduced in Section |2.3j . More 
precisely, if 5 a S(t) is a Gibbs measure for any a, then logpt(o~, rj) is the Hamiltonian of this 
Gibbs measure. In order to prove or disprove Gibbsianness of the measure vS(t), one has to 



study the Hamiltonian ( 3.32 ) for fixed rj. Let us denote by G(H*) the set of Gibbs measures 

Ht(-,rj). From |]i~0|| we have the following: 



associated with the Hamiltonian if* 



Proposition 3.33 For any t > 0: 
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1. If \Q(H*)\ = 1 for all r] E ft, then uS(t) is a Gibbs measure. 

2. For monotone specifications, if \Q{H^)\ > 2, then rj is a bad configuration for vS(t), so 



vS(t) is not a Gibbs measure (by Proposition 2. It ). 



Proof. See [jlCfl . Part 2 is expected to be true without the requirement of monotonicity but 
this has not been proved. I 

A monotone specification arises e.g. when the Hamiltonian of ( |3.32| ) comes from a ferromag- 
netic pair potential and an arbitrary single-site part (possibly an inhomogeneous magnetic 
field). 

In the case of a high-temperature dynamics (T M >> 1), 5 a S(t) converges to \x for any a. 
This implies that for large t we can view the Hamiltonian of ( |3.32| ) as follows: 

H t (a, rj) = H u (a) + H^rj) + o a , v (t), (3.34) 

where o CTiT) (t) means some Hamiltonian with corresponding interaction converging to zero as 
t f oo in B. Therefore, if H u does not have a phase transition, then should not have a 
phase transition either for large t. On the other hand, if H v does have a phase transition, 
then the o<j ir) (t)-term will be important to select one of the phases. In Sections |s| — |o| we will 
come back to this description in more detail. 

The case of independent spin flips corresponds to = 0. 



4 Conservation of Gibbsianness for small times 

Having put the technical machinery in place in Sections [2]-|3], we are now ready to formulate 
and prove our main results in Sections ||-||. 

In this section we prove that for every finite-range spin-flip dynamics starting from a Gibbs 
measure v corresponding to a finite-range interaction the measure uS(t) remains Gibbsian in 
a small interval of time [0, to]- The intuition behind this theorem is that for small times the 
set of sites where a spin flip has occurred consists of "small islands" that are far apart in a 
"sea" of sites where no spin flip has occurred. This means that sites that are far apart have 
more or less disjoint histories. 



Theorem 4.1 Let both the initial measure v and the reversible measure (i be Gibbs measures 
for finite-range interactions XJ V resp. XJn- Then there exists to = to(n,u) > such that vS(t) 
is a Gibbs measure for all < t < to. 

Proof. During the proof we abbreviate = H^' u . We prove that the limit 

^ipm-m (4 . 2) 

E CT (exp[H A (a t )}) v ; 

converges uniformly in t G [0, to] f° r *0 small enough when U^,U^ E B^ r . The to depends on 
both U„ and U^. 

Let us write R u , to denote the range of U v , (see Section p^ ). 
I: R u < oo, R fl = 0. 
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To warm up, we first deal with unbiased independent spin-flip dynamics. For this dynamics 
the distribution of at under P^ coincides with the distribution of af under P^. Therefore we 
can write 

E° (ex P [# A (a?)]) J2 AcA 5[ Al (l-5 t )W-\ A \eM(H AA ^-H)(a)] 



K (exp[£T A (<7t)]) J2acA 4 A '(1 " <^) |AHA| exp([H A - H)(a)} 

CacA (A) exp[(ff^W - 
V Z A cA{A) lAl ^{(H A -H)(a)} 



*x(<r), (4.3) 



where 



* x (a) = exp[(HW -H)(a)] (4.4) 
is a continuous function of a, the sum runs over 

A = {y e A : a t (y) / a (y)}, (4.5) 

while 

<5t = K(?t(x) + <ro(?)) = 1 " e~ 2t . (4.6) 
The notation H , A C A, is defined by 

i^(<r) = ff(<7 A ) (4.7) 

with a" 4 the configuration obtained from a by flipping all the spins in A. 
Suppose first that R u = 1. Then 

h Aub_ h a = h b_ h y A,B : d(A,B) > 1. (4.8) 

For i C A we can decompose A into disjoint nearest-neighbor connected subsets 71,. . . ,7^ 
and thus rewrite (14.3!) as follows: 



E° (exp[ff A (of )]) _ ^«=o h E 7l ,..., 7nC A, 7l n 7j =0nr=i <(7i)^ ^ 



E° (exp[i? A (a t )]) ^=0 & E 71l ... >TnC A, 7i n 7 ,=0 EGU 

with 



*( 7 ) = e[ ll exp[H^ A ^(a)- H^(a)} 



(4.10) 



1^(7) = e[^exp[H^(a) - H(a)} 

and et = 5t/(l — 5t). Note that 10^(7) = uv(7) for all 7 that do not contain x. 
Next, since 

\{H^-H){a)\ <| 7 |C (4.11) 

with 

C = 2 sup sup < 00, (4.12) 

A <r |A| 

we have the estimate 

\wa{j)\ < exp(-a t |7|) with a t = -C + log(l/e t ). (4.13) 
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A similar estimate holds for |u;J(7)|. Since at | oo as t [ 0, it follows that for t small enough we 
can expand the logarithm of both the numerator and the denominator in (^^) in a uniformly 
convergent cluster expansion: 




7i,...,7„CA,7 l ri7j=0 i=l 
n 

^n=0 "" 7 1 ,...,7„cA,7 l ri7j=0 i=l 

By the estimate ( fl.13 ) we have, for t small enough, 



E fl ( r W(r), 

r 



5>(TK(T). 



lim sup 



sup|a(r)«(r)-«; ff (r))| =0 VxG 



(4.14) 



(4.15) 



and hence we obtain uniform convergence of the limit in (4.2). 

The case R v < oo is treated in the same way. We only have to redefine the jiS as the 
.Reconnected decomposition of A. Note that to depends on R u and converges to zero when 
R v | oo. 

II: R u < oo, Ru < oo. 



Next we prove that the limit (4.2) converges uniformly if both interactions U^,U U are finite 
range. For the sake of notational simplicity we first restrict ourselves to the case R v = R^ = 1. 

We abbreviate U = — U u . The idea is that we go back to the independent spin-flip 
dynamics via Girsanov's formula. After that we can again set up a cluster expansion, which 
includes additional factors in the weights due to the dynamics. 

The first step is to rewrite (|4.2|) in terms of the independent spin-flip dynamics: 



E g * (exp[H A (a t )]) 
E a {exp[H A (a t )}) 



(4.16) 



K (exp ( Ey^A Jo log <y, v X s)dN y s + /„*(! - c(y, a*))ds) exp[ff A (a?)] 



K ( ex P [EyeA So lo S c (y' (T s)dNl + / * (1 - c(y, a s ))dsj exp[H A (a t )] 

For a given realization uj of the independent spin-flip process, we say that a site y is u-active 
if the spin at that site has flipped at least once. The set of all w-active sites is denoted by 
J(uj). Let a denote the trajectory that stays fixed at a over the time interval [0, t]. For A C A, 
define 



U x {A,u) 
U 2 {A,uj) 



f \ogc(y,uj s )dN$(uj) + / (1 - c(y,u s ))ds if A = D c 
HA^D C 



U(A,ut) 



and put 
Also define 



U(A,u>) =U 1 (A,lo) +U 2 (A,uj) 
U x (A,u)=U(A,u; x ), 



(4.17) 

(4.18) 
(4.19) 
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where the trajectory u x is defined as 

(u x ) a = (oj s ) x 0<s<t. (4.20) 
With this notation we can rewrite the right-hand side of ( [4.16| ) as 
/ E° (exp (EacaI^(Aw) - U*(A,a)]))' 



(4.21) 



where 



^ E° (exp(£ AcA [W(A,u;)-W(A,<7)])) 

* x (a)=exp i^2[U(A,a)-U(A,a x )]\ (4.22) 

\A3x J 

is a continuous function of a. In order to obtain the uniform convergence of (|4.2D, it suffices 



now to prove the uniform convergence of the expression between brackets in ( 4.21 ). 

As in part I, we decompose the set of u;-active sites into disjoint nearest-neighbor connected 
sets 71, . . . , 7^ and rewrite, using the product character of E^, 

K (exp (£ ACA [W*(A u>) - U X (A, a)})) 



n (exp {^ AcA [U(A,u)-U(A,a)])) 

Xm=0 n! X/7i,...,7„cA,7jn7i=0 Ili=l w a(li 



En=0 n! S7 1 ,...,7„cA,7in7 J =0 Ili=l w a(li 

The cluster weights are now given by 



(4.23) 



w, 



■ ( 7 ) = e^lE° [ 1{J( W ) D 7} exp [ £ o> 7 <7 AVy ) - U(A, a)} ) ) , (4.24) 

. AH7^0 



and an analogous expression for w x after we replace U by Z/P. The factor e*' 7 ' arises from the 
probability 

P° ( J(w) c D A \ U»7i) = e-*l A \ u ^l = e~'l A l JJ e 4 ^l . (4.25) 

i 

Having arrived at this point, we can proceed as in the case of the independent spin- flip 
dynamics. Namely, we estimate the weights w a and prove that 

Wa(l) < e~ athl (4.26) 

with at t oo as t J 0. To obtain this estimate, note that 

K(J(u J )D7)<(l-e- t p. (4.27) 



Then apply to ( [4.24 ) Cauchy-Schwarz, the bounds e < c y < M on the flip rates, and the 
estimate 

C = supsup— - \U(A,a)\ < oo, (4.28) 



A a |A| AnA^0 



to obtain 

Wail) <e Ktbl (l-e^r^ bl f or some K = K(e,M,C). (4.29) 
This clearly implies ( 4.26| ). 



The case R u ,Ru < oo is straightforward after redefining the 7j's. 
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5 Infinite-temperature dynamics 



5.1 Set-up 



In this section we consider the evolution of a Gibbs measure v under a product dynamics, i.e., 
the flip rates c(x, a) depend only on a(x). The associated process {at : t > 0} is a product of 
independent Markov chains on {—!,+!}: 



where T a t x \ is the Markov chain on {—1, +1} with generator 

L x (p(a) = c(x,a)[tp(-a) - (p(a)]. 



(5.1) 



(5.2) 



Let us denote by pf(a,{3) the probability for this Markov chain to go from a to (5 in time t. 
The Hamiltonian ( 3.32j ) of the joint distribution of (ctq, o~t) is then given by 



(5.3) 



This equation can be rewritten as 

H t (a, V ) = H u (a) -^^(t)a(x) - ^^(^(x) -J^^M^ 



(5.4) 



with 



hf(t) 

W) 

hf 2 (t) 



X 








Pt(-, 




pU+, 


+K(-,+) 


p* (+. 


-)pf(->-) 


pU+, 






-)pf(-i+) 



(5.5) 



The fields hf resp. /if tend to pull a resp. r] in their direction, while hf 2 is a coupling between 
ex and r/ that tends to align them. Indeed, note that hf 2 (t) is positive because 



Pt(+,+)Pt(-,-)-pU+,-)pU-,+) = det(e W (tL x )) > 0. 



(5.6) 



In what follows we will consider the case where the single-site generators L x are independent 
of x and are given by 



2 



-1 + e 1-e 
1+e -1-e 



for some < e < 1. 



(5.7) 



For e > this means independent spin flips favoring plus spins, for e = it means independent 
unbiased spin flips. The invariant measure of the single-site Markov chain is (^(+), v{— )) = 
5(l + e, 1 — e). The relevant parameter in what follows is 



1 



1 + e' 



(5.8) 
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In terms of this parameter the fields in (|5.5|) become 



h 2 (t) 



1 



\ogS + hi(t) 



M*) = i lo s (1 _ e -t )2 

In particular, for 5 = 1 we get /ii(t) = /^(i) = and 

1 1 + e~* 



(5.9) 



(5.10) 



5.2 



1 << T„ < oo, = oo 



Theorem 5.11 -Let v be a high- or infinite-temperature Gibbs measure, i.e., its interaction 
U v satisfies ^2.2$) . Let S(t) be the semigroup of an arbitrary infinite-temperature dynamics. 
Then uS{t) is a Gibbs measure for all t > 0. 



Proof. The joint distribution of (co,0t) is Gibbs with Hamiltonian (recall ( |3.32| ) and (|5.4|)) 
H t (a, r,) = H v {a) + £[M*) + h l2 (t) V {x))a{x) + h 2 (t) ^ Vx . (5.12) 

X X 

For fixed r/, the last term is constant in a and can therefore be forgotten. Since Ht(-,r/) 



differs from H u (-) only in the single-site interaction, Ht(-,rj) satisfies ( 2.2S| ) if and only if 
#„(•) satisfies ( p8|) . Hence \g(H t (-,rj)\ = 1 for any 77, and we conclude from Proposition |3.33| 
that vS(t) is Gibbsian. I 



Theorem 5.11 should not come as a surprise: the infinite-temperature dynamics act as a 



single-site Kadanoff transformation and in the Dobrushin uniqueness regime such renormalized 
measures stay Gibbsian \14\, JTSfl , ||. 



5.3 < T v « 1, T M = 00, 6 = 1 

For the initial measure we choose the low-temperature plus-phase of the ci-dimensional Ising 
model, v = vphi i- e -) the Hamiltonian H u is specified to be 



H u (a) = -/3j2 < X )<V) - h Y,< a 



(5.13) 



<x,y> 



where ^2 <x y> denotes the sum over nearest-neighbor pairs, and j3 » (3 C with f3 c the critical 
inverse temperature. The dynamics has generator 



Lf = ^V x f, 



(5.14) 



corresponding to the case 5 = 1. The joint measure has Hamiltonian as in ( glD , with 
h^t) = h 2 (t) = and h 12 (t) = h t : 

<x,y> x x 

The "dynamical field" is given by ht = — (1/2) log[tanh(i/2)]. 
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Theorem 5.16 For j3 » j3 c : 

1. There exists a to = to{(3, h) such that vp^S{t) is a Gibbs measure for all < t < to. 



2. If h > 0, then there exists a t\ 
t>h. 



ti((3,h) such that ushS(t) is a Gibbs measure for all 



3. If h = 0, then there exists a t 2 = i2 (/?) such that v/3fiS(t) is not a Gibbs measure for all 
t > t 2 . 

4- For d > 3, if h < h{(3) small enough, then there exist t% = t^{(5,h) and t^ = t^{(5,h) 
such that vp^S^) is not a Gibbs measure for all t% < t < 

Proof. The proof uses (|5.15| ). 

1. For small t the dynamical field ht is large and, for given rj, forces a in the direction of r\. 
Rewrite the joint Hamiltonian in (5.15) as 



JL 



<x,y> 



h t H t (a, rj). 



(5.17) 



For < t < to small enough, Ht has the unique ground state n and so, for A > Ao large 
enough, XH t satisfies (|2.28|) (see [^31 , example 2, p. 147). Therefore, for < t < t\ such 
that \fh~t > Ao, Ht(-,7]) has a unique Gibbs measure for any rj. Hence, z/S^i) is Gibbs by 
Proposition |3.33| (1). 

2. For large t the dynamical field ht is small and cannot cancel the effect of the external field 
h > 0. Rewrite the joint Hamiltonian as 

= y/pH t (tT,ri). (5.18) 

For t>t\ large enough (independently of (3), H t (-,r]) has the unique ground state a = h/\h\. 
Hence, for (5 large enough, y/j3H t (-,r]) has a unique Gibbs measure by ( 2.28 ) (again, see [|b|, 
example 2, p. 147). Hence, uS(t) is Gibbs by Proposition 3.33| (1), 

3. This fact is a consequence of the results in [9], section 4.3.4, for the single-site Kadanoff 
transformation. Since the joint Hamiltonian in ( 5.15| ) is ferromagnetic, it suffices to show that 
there is a special configuration rj spec such that \G (H (■ , rj spec )\ > 2. We choose rj spec to be the 
alternating configuration. For t > t 2 large enough, Ht(-,rj spec ) has two ground states, and 
by an application of Pirogov-Sinai theory (see (9j Appendix B), it follows that, for (3 large 
enough, \G(Ht(-,rj spec )\ > 2. Therefore n spec is a bad configuration for vS(t), implying that 
vS{t) is not Gibbs by Proposition 3.33] (2). 

4. In this case we rewrite the Hamiltonian in ( 5.15[) as 

H t (a,r)) = -I3 Y o-(x)a(y) ~Yi h + h Mx)]a(x). 



(5.19) 



<x,y> 
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For "intermediate" t we have that h and ht are of the same order. As explained in M 
section 4.3.6, we can find a bad configuration rj spec such that the term ^ x htT](x)a(x) in the 
Hamiltonian "compensates" the effect of the homogeneous-field term ^ ho~(x) and for which 
Ht(- , Tjspec) has two ground states which are predominatly plus and minus. Since the proof of 
existence of r] spec requires analysis of the random field Ising model, we have to restrict to the 
case d > 3 (unlike the previous case rj spec is not constructed, but chosen from a measure one 
set). Then for (3 large enough, by a Pirogov-Sinai argument (see appendix B, Theorem B 31 
of §§) \g(H t (-,r] spec )\ > 2, implying that uS(t) is not Gibbs by Proposition p r 33|(2). I 

Remark: 

From the estimate (B89) in Q, Appendix B, we can conclude the following: 

1. to((3, h) — > as j3 — ► oo, and to(f3, h) — ► oo as h — ► oo. 

2. t 2 (/3) -> as -> oo. 

3. tz((3, h) — ► as j3 — > oo. 

5.4 < T v « 1, Tf, = oo, 5 < 1 

Let us now consider a biased dynamics. At first sight one might expect this case to be 
analogous to the case of an unbiased dynamics with an initial measure having h > 0. However, 
this intuition is false. 



Theorem 5.20 The same results as in Theorem 5. It hold, but with the ti's also depending 
on 5. For item 4 we need the restrictions d > 3 and \h + \ log <5| small enough. 



Proof. The last term in ( p.4[ ) being irrelevant, we can drop it and study the Hamiltonian 

H t (a, V ) = -(3 o-(x)a(y)-Y,^)[(h + h 1 (t)) + h 12 (t)7 1 (x)}. (5.21) 



<x,y> 



This Hamiltonian is of the same form as ( 5.15| ), but with h becoming i-dependent. We have 
lim^oo h\(t) = and lim^oo h\2 if) = with 

lim Ma = i±| >1 , (5.22) 

itoo hi(t) 1-6 V ' 

so that, in the regime where (3 » {3 C , h = 0, t >> 1, we find that the effect of h^ft) 
dominates. Hence we can find a special configuration that compensates the effect of the field 
h\(t) and for which the Hamiltonian ( 5.21| ) has two ground states, implying that uS(t) Q. 



Similarly, when h > we can find t intermediate such that 'Yl ix {h+hi{t))a{x) is "compensated" 
by J2 h i2(t)o-(x)rj(x). I 

Remark: 

Note that if T v = 0, T^ = oo, then uS(t) is a product measure for all t > and hence is Gibbs. 
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6 High-temperature dynamics 



6.1 Set-up 

In this section we generalize our results in Section |E| for the infinite-temperature dynamics to 
the case of a high-temperature dynamics. The key technical tool is a cluster expansion that 
allows us to obtain Gibbsianness of the joint distribution of (o"Q,cr t ) with a Hamiltonian of 
the form ( 3.32; ). The main difficulty is to give meaning to the term logpt(<r, 77), i.e., to obtain 



Gibbsianness of the measure 5 a S(t) for any a. In the whole of this section we will assume 
that the rates c(x, a) satisfy the conditions in Section |2.2j and, in addition, 



c(x, a) = 1 + e(x, a) (6-1) 

with 

snp atX \e(x,a)\ = 5 «l , g ^ 

e(x, a) = e(x, —a). 

The latter corresponds to a high-temperature unbiased dynamics, i.e., a small unbiased per- 
turbation of the unbiased independent spin-flip process. For the initial measure we consider 
two cases: 

1. A high- or infinite-temperature Gibbs measure v. In that case we will find that uS(t) 
stays Gibbsian for all t > 0. 

2. The plus-phase of the low-non-zero-temperature d-dimensional Ising model, v$ h) corre- 
sponding to the Hamiltonian in Q5.13j ). In that case we will find the same transitions as 
for the infinite-temperature dynamics. 

6.2 Representation of the joint Hamiltonian 

In this section we formulate the main result of the space-time cluster expansion in p8[ and 



[25|. We indicate the line of proof of this result, and refer the reader to [25] for the complete 



details. 

Theorem 6.3 Let v be a Gibbs measure with Hamiltonian H u , and let the dynamics be gov- 



erned by rates satisfying (6.1- 6. i ). Then the joint distribution of (o"o,0t), when ctq is dis- 



tributed according to v , is a Gibbs measure with Hamiltonian 

H t (a, V )=H l/ (a)+H t dyn (a,r ] ). (6.4) 

The Hamiltonian H^ yn (a,r]) corresponds to an interaction U dyn (A,a,r]), A £ S, that has the 
following properties: 

1. The interaction splits into two terms 

U\yn = ^0 + Ul (6.5) 
where Uq is the single-site potential corresponding to the Kadanoff transformation: 

U&({x},a,ri) = -ilog[tanh(t/2)]a(x)7?(x) x£Z d , 
U^(A,a,rj) = if \A\ ± 1. 
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2. The term U\ = Ug(A,o~,<n) decays exponentially in the diameter of A, i.e., there exists 
a(6) > such that 

supsup Vsupe a(5)diam(A) |C/|(A )f 7,?7)| < oo. (6.7) 
and a(S) | oo as 5 j 0. 

3. The potential converges exponentially fast to the potential of the high-temperature 
reversible Gibbs measure: 

limsup V sup e a(5)diam(A) | ^(Aa,!?) - U^rj)] = 0. (6.8) 

4- The term Ut is a perturbation of the term Uq, i.e., 

v sup x J2 A Bx sn Pay, v \Us(A,a,r]) - U^A,?' ,n)\ _ 

lim SUP : ; — ; — — U. I O.i) ) 

6io t >5 log[tanh(t/2)] 

Remarks: 

1. Equation (|6.6|) corresponds to the infinite-temperature dynamics (i.e., c = 1). 



2. Equation (|6.9| ) expresses that the potential as a function of the rates c is continuous at 
the point c = 1, and that the Kadanoff term is dominant for 5 « 1. 



Main steps in the proof of Theorem 6.3 in [25]: 



Discretization: The semigroup S(t) can be approximated in a strong sense by discrete- 
time probabilistic cellular automata with transition operators of the form P n (o~' \o~) = 
H x p n(o-'(x)\a), where 



Pn(o-'(x)\a) 



1 c(x, a) 

n 



Space-time cluster expansion for fixed discretization n: For n fixed the quantity 

= log ^^Sr (6.11) 

is defined by the convergent cluster expansion 

K(°,V)= E (6-12) 
TBx,reC 

where C is an appropriate set of clusters on Z rf+1 . 

Uniformity in the discretization n: The functions converge uniformly as n f oo 
to a continuous function ^ x (which defines a continuous version of dfj, x /djj). This is 
shown in two steps: 
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1. Uniform boundedness: 



sup sup sup I ^ (a, rj) \ < oo. 

n x a,T) 



(6.13) 



2. Uniform continuity: 



lim supsupsup|1^(<7ACAe,?7A£A<0 - *£0,r/)| =0 Vo", 77 G $7. (6.14) 



Equations ( 6.13| ) and (6.14) imply that ^ as a function of n contains a uniformly 
convergent subsequence. The limiting fy x is independent of the subsequence, since it is 
a continuous version of dfi x /d(i. 

6.3 1 << T v < 00, 1 << Tu. < 00 



Given the result of Theorem 6.3, the case of a high- or infinite-temperature initial measure is 



dealt with via Dobrushin's uniqueness criterion (recall Theorem 5.11 in Section |5.2| ). 



Theorem 6.15 Let v be a high-temperature Gibbs measure, i.e., its interaction U v satisfies 
\2.2t ). Let the rates satisfy \6.1-6.i). Then, for 5 small enough, vS(t) is a Gibbs measure 
for all t > 0. 



Proof. For fixed rj, the Hamiltonian Ht(-,r)) of d6.4| ) corresponds to an interaction U^' . By 
(|6.7D and ( |6.9| ), this interaction satisfies 



limsupsup^(|A| - l)sup|l^V) " Uf{a') 



A3x 



-£(\A\-l)\U u (a)-U u (o J )\<2. 



(6.16) 



Therefore, for 5 small enough, ( p. 28 ) is satisfied for the interaction U^' S for all t > and all 
rj. Hence \Q(H t {-,rj))\ = 1, and we conclude from Proposition 3.33| (1) that vS{t) G Q. I 



6.4 < T v « 1, 1 « Tu < 00. 

We consider as the initial measure the plus-phase of the low-temperature Ising model vphi 



introduced in Section 5.3. The joint distribution of (00,04) has the Hamiltonian 
H t (a, V ) = -P ^(x)a(y)-/ l ^a(x)-^log[tanh(t/2)]^a(x)r / (x)+^(0,r / ), (6.17) 

<x,y> x x 

where H$ corresponds to the interaction E/J introduced in ( |6 . 5[ ) . The following is the ana- 
logue of Theorem 5.16| 



Theorem 6.18 For (3 » f3 c and < 5 « 1: 

1. There exists to = to(/?, h, 5) such that ug^Sif) is a Gibbs measure for all < t < to- 
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2. If h > 0, then there exists t\ = t\(f3,h,§) such that i/QhS(t) is a Gibbs measure for all 
t > ti. 

3. If h = 0, then there exists ti = t2((3,5) such that upflS{t) is not a Gibbs measure for all 
t > t 2 . 

4- For d > 3, if < h < h{j3) and < 6 < 5((3,h), then there exists t${(3,h,S),ti{(i,h,$) 
such that vafrSit) is not a Gibbs measure for all t £ [£3, £4]. 

Proof. 

1. This a consequence of Theorem 4.1. 



2. This is proved in exactly the same way as the corresponding point in Theorem 5.16 



3. Here we cannot rely on monotonicity as was the case in Theorem 5.16. It is therefore 
not sufficient to show that for the fully alternating configuration rj a , the Hamiltonian 
H(-,rj a ) exhibits a phase transition. We have to show the following slightly stronger fact: 
if m\{do~) is any Gibbs measure corresponding to the interaction H(-, ?7^+a<0, then 

J ml(da)a(0) > 7 > 0. (6.19) 

This proof of this fact relies on Pirogov-Sinai theory for the Hamiltonian Hf(-, r^+A )- 
The first step is to prove that the all-plus-configuration is the unique ground state of this 
Hamiltonian. Since the Ising Hamiltonian satisfies the Peierls condition, we conclude 
from ||] Proposition B.24 that the set of ground states of H t (- , n^+jj is a subset of 
{+, — }. If we drop the term H^(-, ?7^+a c ) (he., if 5 = 0), then the remaining Hamilto- 
nian has as the unique ground state the all-plus-configuration and satisfies the Peierls 
condition. Therefore, for 5 small enough, we conclude from Q Proposition B.24 that 
Ht('i ^a+A 2 ) has the all-plus-configuration as the only possible ground state. From ( |6.17| ) 
it is easy to verify that the all-plus-configuration is actually a ground state for 5 small 
enough. In order to conclude that for (3 large enough, the unique phase of i?t(-, ?7a+A c ) 
is a weak perturbation of the all plus configuration (uniformly in A), we can rely on 
the theory developed in Q, or || which allows exponentially decaying perturbations of 
a finite range interaction satisfying the Peierls condition (see e.g. equations (1.3), (2. 2) 
of [^]). Similarly, Ht(- , r]A~ A c ) has a unique phase which is a weak perturbation of the 
all minus configuration. This is sufficient to conclude that no version of the conditional 
probabilities is continuous at r] a , see the discussion [|[ p. 980-981. 

4. We can use the same argument as developed in ||, section 4.3.6, introducing a random 
perturbation of the alternating configuration to "compensate the uniform magnetic field" 
(since this requires analysis of the random field Ising model, we have the restriction 
d > 3). The only complication is the extra term in the Hamiltonian arising from 
5^0. This requires Pirogov Sinai theory for the interaction Ht(o~,r]), where rj = if 
is a random modification of the fully alternating configuration obtained by flipping the 
spins in the alternating configuration with probability e/2/3 for a flip from + to — . Since 
the couplings between n and a are not finite range, we cannot apply directly Theorem 
B31 of H for the random Hamiltonian Ht(o-,r] e ). However, as the interaction decays 
exponentially fast and Pirogov-Sinai analyses do not distinguish between finite range 
and exponentially decaying interactions, similar arguments as those developed in EH] 
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still work in our case and yield the analogue of Theorem B31 of Ml. However we have 
not written out the details. 



Remark: 

A result related to Theorem |6.15 was obtained in |2£| . Although the abstract of that paper is 
formulated in a somewhat ambiguous manner, its results apply only to initial measures which 
are product measures (in particular Dirac measures) . In particular this includes the case 
T v = and 1 << T M < oo. The results of |^8| (or j25fl ) then imply that the measure is Gibbs 
for all t > 0. This seems surprising, because t2(P,5) j as (3 j oo. It is therefore better for 
the intuition to imagine a Dirac-measure as a product measure than to view it as a limit of 
low-temperature measures. 



7 Discussion 

7.1 Dynamical interpretation 

In the case of renormalization-group pathologies, the interpretation of non-Gibbsianness is 
usually the presence of a hidden phase transition in the original system conditioned on the 
image spins (the constrained system). In the context of the present paper, we would like to 
view the phenomenon of transition from Gibbs to non-Gibbs as a change in the choice of most 
probable history of an improbable configuration at time t > 0. 

To that end, let us consider the case of the low-temperature plus-phase of the Ising model 
in zero magnetic field (/? >> f3 c , h = 0) with an unbiased (5 = 1) infinite-temperature 
dynamics. Consider the spin at the origin at time t conditioned on a neutral (say alternating) 
configuration in a sufficiently large annulus A around it. For small times the occurrence of 
such an improbable configuration indicates that with overwhelming probability a configuration 
very similar was present already at time 0. As the initial measure is an Ising Gibbs measure, 
the distribution at time of the spin at the origin is determined by its local environment only 
and does not depend on what happens outside the annulus A. As all spins flip independently, 
no such dependence can appear within small times. 

However, after a sufficient amount of time (larger than the transition time t2), if the 
same improbable configuration is observed, then it has much more chance of being recently 
created (due to atypical fluctuations in the spin-flip processes) than of being the survivor of 
an initial state of affairs. Indeed, to have been there at time is improbable, but to have 
survived for a large time is even more improbable. Suppose now that outside the annulus A we 
observe an enormous annulus T in which the magnetization is more negative than —m*(t)/2, 
where m*(t) is the value of the evolved magnetization (which starts from m*(0) and decays 
exponentially fast to zero). Because a large droplet of the minus-phase shrinks only at finite 
speed and typically carries a magnetization characteristic of the evolved minus-phase, with 
large probability there was an enormous droplet of the minus phase (even a bit larger than T) 
at time 0, which the spin at the origin remembers. Indeed, the probability of this happening 
is governed by the size of the surface of T. In contrast, the probability of a large negatively 
magnetized droplet, arising through a large fluctuation in the spin-flip process starting from 
a typical plus-phase configuration, is governed by the volume of T. Therefore, this second 
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scenario can safely be forgotten. Although for any size of the initial droplet of the minus- 
phase there is a time after which it has shrunk away, for each fixed time t we can choose an 
initial droplet size such that at time t it has shrunk no more than to size T. Since we want the 
shrinkage until time t to be negligible with respect to the linear size of V, we need to choose 
r larger when t is larger. 

Thus, the transition reflects a changeover between two improbable histories for seeing an 
improbable (alternating) annulus configuration. It can be viewed as a kind of large deviation 
phenomenon for a time-inhomogeneous system. One could alternatively describe it by saying 
that for small times a large alternating droplet must have occurred at time 0, while after the 
transition time £2 a large alternating droplet must have been created by the random spin- flips: 
a "nature to nurture " transition |35j . The mathematical analysis of this interpretation would 
rely on finding the (constrained) minimum of an entropy function on the space of trajectories. 
Alternatively, one could try to study the large deviation rate function for the magnetization 
of the measure at time conditioned on an alternating configuration at time t. This rate 
function should exhibit a unique minimum for < t < £2 and two minima for t > £2- 



7.2 Large deviations 

A measure can be non-Gibbsian for different reasons (see ||, section 4.5.5) One of the possi- 
bilities is having "wrong large deviations", i.e., the probability 



uS(t) 



r x f(a) 



a 



(7.1) 



for fixed t and a / J S(t)fdv does not decay exponentially in |A|, i.e., not as exp[— |A|/j(q) + 
o(|A|)], or equivalently, there exists a function / e C, f > 0, / / such that 



bm n-log / uS(t)(da)exp 
ATZ d |A| J 



r x f(cr) 



(7.2) 



An example where this phenomenon of "wrong large deviations" occurs is the stationary 
measure of the voter model (see e.g. p2]]). However, it does not occur in our setting. Namely, 
if the scale of the large deviations of the random measure L\ = X^ceA 5 Tx(T under v is the 
volume |A|, then the same holds under uS(t) for any t > 0. Indeed, by Jensen's inequality 
and by the translation invariance of the dynamics we have, for / E C, f > 0, / 7^ 0, 



1 f 

lim sup -j— -r log / uS(t)(da) exp 
AfZ d l A l J 

> lim sup — -7 log / v (da) exp 
ATZ d \ A \ J 

S(t)fdn-h(jj,\u) 



.xeA 



sup 



> 



(7.3) 



with h(-\-) denoting relative entropy density. The equality follows from the volume-scale large 
deviations of v, and the last inequality follows from the fact that S(t)f € C(O), S(t)f > 0, 
S(t)f / imply / S(t)fdu > 0. 
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7.3 Reversibility 



Throughout the whole paper, we have assumed the stationary measure [i to be reversible. 
However, this is a condition that only serves to make formulas nicer. It is not at all a 
necessary condition: if we consider any high-temperature spin-flip dynamics, then we know 
that the stationary measure \i is a high-temperature Gibbs-measure. Equation ( [3.2; ) can be 
rewritten in the general situation: we have to replace S\(t) in the right-hand side by S^(t), 
where S*(t) is the semigroup corresponding to the rates of the reversed process, i.e., the rates 

c*(x,a)=c(x,a x )^. (7.4) 

In all the formulas of Section 2, we then have to replace E CT by E* , referring to expectation in 
the process with semigroup S*(t). 



7.4 Open problems 

1. Infinite-range interactions. How much can we save when relaxing the condition that 
the interactions be finite-range? 

2. Trajectory of the interaction. In the regime 1 << T v < oo, 1 << < oo, what 
can we say about the trajectory t \— ► lit? It is not hard to prove that it is analytic 
in Bu and converges to U^. But can we say something about the rate of convergence? 
Note that we can view the curve {U Ut : t > 0} as a continuous trajectory in the space 
B, interpolating between XJ V and Un, which implies that Q contains an arc-connected 
subset. Other topological characteristics of Q are discussed in |J, section 4.5.6. 

3. Uniqueness of the transitions. Even in the case = oo we have not proved that 
the transition from Gibbs to non-Gibbs is unique e.g. that to(/3, 0) = ^(Z?) in Theorem 
5.16] , However, we expect that when h = the alternating configuration is "the worst 



configuration", i.e., the transition is sharp and occurs at the first time at which the 
alternating configuration is bad. 

Estimates for the transition times. Can we find good estimates for the ti's as a 
function of e.g. the temperatures, the magnetic fields and the ranges of the interaction 
in v and \x. 

Weak Gibbsianness. In the regimes where vS{t) is not a Gibbs measure we expect that 
we can still define a vS(t)-&.s. converging interaction U% for which vS(t) is a "weakly 
Gibbsian measure" (see |J, ^7j). This interaction Ut can e.g. be constructed along 



similar lines as are followed in the proof of Kozlov's theorem (see [21|,[pq]) and its 
convergence is to be controlled by the decay of "quenched correlations", i.e., the decay 
of correlations in the measure at time conditioned on having a fixed configuration ij at 
time t. These correlations are expected to decay exponentially for uS(t)-a,.e. r/, which 
would lead to vS(t)-&.s. convergence of the Kozlov-potential. 

Low-temperature dynamics. The main problem of analyzing the regime < T M << 
1 for large t is the impossibility of a perturbative representation of — log pt(cr, rf). If we 
still continue to work with the picture of the joint Hamiltonian in ( |3.32 ), then the term 



— logger, rf) will not converge to a cr-independent Hamiltonian as t f oo. Therefore we 
cannot argue that for large t the Gibbsianness of the measure uS(t) depends only on the 
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presence or absence of a phase transition in the Hamiltonian H u of the initial measure 
v. The dynamical part of the joint Hamiltonian can induce a phase transition. The 
regime < « 1 is very delicate and there is no reason to expect a robust result for 
general models. Metastability will enter. 

7. Zero-temperature dynamics. What happens when = 0? In this case there is only 
nature, no nurture. We therefore expect the behavior to be different from < T M << 1. 
Trapping phenomena will enter. 

8. Other dynamics. Do similar phenomena occur under spin-exchange dynamics, like 
Kawasaki dynamics ? In particular, how do conservation laws influence the picture (see 

a 0, 0)? 
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